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ABS TRACT 

) — 

The  model  developed  in  this  report  is  an  extension  and 

reformulation  of  a model  called  the  Coherence  model  for 

guiding  EEO  (Equal  Employment  Opportunity)  planning  at  the 

micro-level  in  the  U.S.  Navy's  civilian  workforce  developed 

by  Charnes,  Cooper,  Lewis  and  Niehaus.  This  model  Is  called 

cue  Goal-Arc  model. 

Like  its  predecessors,  the  Goal-Arc  model  utilizes  a 
goal  programming  approach  with  embedded  Markoff  processes. 

As  in  the  Coherence  model,  piecewise  linear  goal  functionals 
with  "artifact  goals"  are  used  to  approximate  the  transition 
relations  of  the  Markoff  process.  The  Goal-Arc  model,  however, 
carries  this  to  another  stage  of  development.  Analytical  as 
well  as  network  formations  and  interpretations  are  provided 
in  the  following  article.  A numerical  example  with  related 
interpretations  for  EEO  planning  is  also  provided. 


Introduction 


The  Multi-Level  Coherence  Model  for  Equal  Employment 
Opportunities  (EEO)  planning  of  Charnes,  Cooper,  Lewis  and  Niehaus 
(see  [5])  was  developed  in  a dyadic  format.  It  was  less  general  in  its 
development,  however,  than  might  be  required  for  some  cases.  For  many 
applications,  recourse  is  needed  to  large  scale  highly  efficient 
network  codes  such  as  PNET- ^ (which  the  U.S.  Navy  Office  of  Civilian 
Personnel  (OCP)  has)  which  can  readily  handle  multiple  arcs  between 
nodes  as  well  as  lower  and  upper  bounds  on  arc  flows.  The  new,  more 
general,  nonlinear  goal-arc  network  model  discussed  in  this  paper  was 
therefore  developed  to  exploit  these  and  other  possibilities. 


A 


I 


— See  [7]  and  [8].  Other  codes,  such  as  GNET,  are  discussed  in  [lj  and 
[2] . See  also  [9] . 
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In  pa9t  research  papers  (Cas9,  Charnes,  Cooper,  Niehaus  [3]) 
goal  programming  models  of  distribution  (or  assignment)  type  have  been 
reformed  into  equivalent  models  of  distribution  type.  The  "MEEO"  — 
Multi-Level  EEO  — model  was  the  first  to  approximate  Markoff  transition 
constraints  by  "goal  artifacts”  whicn  replaced  the  constraints  by  goals 
with  convex  goal  functionals  on  certain  dyadic  cell  elements.  Here  we 
use  the  analogous  device  for  networks:  the  cells  with  nonlinear  goal 
functions  are  replaced  by  arcs  with  goal  functionals  which  we  shall  now 
call  "goal-arcs."  The  network  format  with  "transhipment"  nodes  which  we 
shall  now  Introduce  allows  us  to  simplify  by  dispensing  with  the 
transshipment  elaboration  and  the  extra  rows  and  columns  this  required 
in  the  dyadic  format.  Finally,  to  fit  the  data  format  of  the  PNET  code 
a "supersource"  and  "supersink"  is  introduced  which  connect  to  arcs  whose 
bounds  replace  the  influxes  and  effuxes  in  other  model  elaborations. 

In  summary,  this  is  part  of  a continuing  evolution 
involving  an  interplay  between  practically  oriented  implementations 
and  research  which  started  at  the  U.S.  Navy,  as  in  [6],  with  a 
combined  goal-programming  Markoff  process  model  for  joining  EEO 
with  civilian  manpower  planning  in  terms  of  targeted  goal  for  each 
of  them.  This  was  followed  by  the  Coherence  («  MEEO)  model  which, 
aa  already  noted,  was  the  first  to  approximate  Markoff  transition 
constraints  with  "goal  artifacts"  in  a suitable  dyadic  formulation. 

See  [5].  Now  we  replace  the  latter  with  a network  model  with 
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i transshipment  nodes  that  make  it  possible  to  obtain  access  to 

I 

currently  available  ultra  high  speed  computer  codes.  This,  in  turn, 
should  make  it  possible  to  provide  interactive  computer  capabilities, 
if  desired,  by  means  of  which  manpower-EEO  planning  can  be  directly 
integrated  into  management  decision  making  instead  of  being  confined 
for  separate  processing  by  customary  personnel  department 
specializations. 

The  way  in  which  this  network  formulation  is  achieved  and 
some  of  the  uses  to  which  it  can  (and  will)  be  put  are  described  in 
the  sections  that  follow.  First  we  shall  describe  some  of  the  node-arc 
conventions  we  employ.  Then  we  shall  provide  an  analytic  characterization. 
The  network  representation  developed  from  this  analytic  model  will  then 
be  depicted  and  the  goal  arc  decompositions  described  to  show  how  the 
convex  functional  elements  are  accommodated. 

This  will  be  followed  by  a numerical  example  which  will  be 
similarly  developed  and  interpreted.  The  resulting  solutions  will 
be  portrayed  in  the  form  of  reports  for  possible  managerial  use 
that  will  help  to  point  up  some  of  these  possibilities  via  the 
prototype  (toy)  example  we  shall  be  employing.  This  will  be 
followed  by  a Summary  and  Conclusion  section  that  will  also  suggest 
some  possibilities  for  further  research. 
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The  Goal-Arc  Model 

We  shall  describe  our  transfers  in  terms  of  flows  on  several 
types  of  arcs  between  several  types  of  nodes: 

(1)  To  each  job  in  each  period  we  assign  two  nodes,  an 
"antecedent"  and  a "consequent."  We  also  designate  as  "job"  nodes 
those  corresponding  to  outside  sources  for  recruitment  (antecedents) 
and  outside  involuntary  retirements  (consequents) . We  also  designate 
"Job  node8"for  normal  organizational  attrition  (consequents).  We 

designate  the  class  of  antecedent  "job"  nodes  for  period  t as  J (t);  the 
class  of  consequent  "job"  nodes  by  J+(t).  (t)  is  the  iC^  job 

antecedent  node;  Jj  + (t)  is  the  j1'*1  job  consequent  node. 

(2)  For  each  proper  (real)  job  between  two  periods  we 
designate  a "valve"  node  to  receive  the  goal  arc  flow  from  the 
consequent  node  of  the  immediate  past  period  and  to  transmit  an 
upper  and  lower  bounded  flow  to  the  next  period  antecedent  node.  We 
let  V^(t)  denote  the  valve  node  for  job  i between  periods  t-1  and  t. 

(3)  A supersource  node,  S , and  a supersink  node,  S . , , 

o m-1 

are  added  for  PNET  code  purposes.  The  supersink  node  is  connected 
back  to  the  supersource  node.  Thereby  every  node  becomes  a tran- 
shipment node. 


i 
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The  flow  on  every  arc  is  unidirectional.  The  arcs  may  be 
"goal"  arcs  (with  a nonlinear  goal  functional)  involving  multiple 
arcs  between  the  same  two  nodes,  or  they  may  be  simple  arcs.  Every 
simple  arc  (or  individual  arc  of  multiple  arcs)  may  have  an  upper  and 

a lower  bound  on  its  flow. 

ic  •>  4* 

Let  x^(t)  denote  the  flow  from  node  (t)  to  node  (t) 

on  the  kC^  individual  arc  of  a multiple  "goal  arc."  The  corresponding 

k k 

lower  and  upper  bounds  are  L^(t)  and  U^j(t). 

Let  xQi  denote  the  flow  from  the  supersource  to  (1).  Let 

x,  , denote  the  flow  from  J.+(n)  to  the  supersink.  Let  x - denote 
i n+1  i n+JL o 

the  flow  from  the  supersink  to  the  supersource. 

Jq 

Let  y^(t)  denote  the  flow  on  arc  k of  the  goal-arc  between 
Ji+(t-l)  and  V^(t) . The  corresponding  upper  and  lower  bounds  are 

tr  If  — 

L^(t)  and  U^(t).  Let  y^(t)  denote  the  flow  on  the  "valve  arc 
between  V^(t)  and  (t). 

The  network  node  conditions  may  now  be  written  explicitly: 


(1)  for  supersource 


x . . - lx 

n+1  o oi 

i£j"(l) 


(2)  for  J (1) 


t.  - I . E x*(l) 
oi  jEJ+(l)  k ij 


(3)  for  JT(1) 


k ieJ  (1) 


«u(1> 


o,  J i J. 


I 
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Where  jQ  is  the  "outside"  node, 


Xoj  + 1 - xij  (1)  " < 

Jo  k ieJ  (1)  1Jo  r Jo 


x^,  (1)  - Zy,  (1)  - 0 


Note  that  there  is  never  flow  from  the  "outside"  node  J~  (1)  to  the  natural 

+ J° 
attrition  node  (t).  We  also  have 

o 

(4)  for  V1(t) 


E y,(t)  - yi(t) 

k 1 1 


(5)  for  (t)  , t>l 


y1(t)  - Z Z x *(t) 

k j 


(6)  for  Jj  (t)  , f>l 


£ £ _ xi?  - E y^(t) 

k ieJ  (t)  1J  r J 


(7)  for  supersink  S 


£ y,  (t)  + Z x . - x - 0 

t *o  ieJ  (n)  1 n 1 0+10 


! 
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An  illustration  of  the  Goal-Arc  Model  is  given  i.i  Figure  1 for  n 

time  periods  and  nrt-2  job  categories.  Sq  is  the  supersource  node  introduced 

on  the  left  and  is  the  supersink  node  introduced  on  the  right.  In 

the  diagram  the  antecedents  and  the  consequents  of  the  outside  node  are 

represented  by  J^Ct)  = J “(t),  J^_  (t)  = J +(t),  and  J^2(t)  ” ^(t). 

^ o ^ o o 


FIGURE  1:  GOAL  ARC  NETWORK  MODEL 
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Some  of  Che  arcs  represent  natural  flows  and  some  may  be 
, goal  arcs.  Recall  that  the  purpose  of  each  of  the  goal  arcs  Is  to 

! represent  a nonlinear  goal  functional  element.  To  represent  these 

piecewise  linear  (nonlinear)  goal  elements  we  can  replace  each  goal 
arc  by  multiple  capacitated  arcs  between  the  same  two  nodes.— ^ 

An  Illustration  is  supplied  in  Figure  2.  The  arc  G between 

1 

nodes  and  Is  a goal  arc.  This  is  indicated  by  the  symbol 
which  we  have  omitted  from  these  links  in  Figure  1 to  avoid  further 
cluttering  of  the  diagram. 

The  lower  portion  of  Figure  2 shows  the  decomposition.  The 

k k k k 

flow  z on  G is  broken  up  into  flows  z on  G where  £z  =z.  Each  z is 

k 

k 

a bounded  variable.  Further  we  let  c be  the  slope  assigned  for  the 
flow  z . Thus, the  decomposition  of  the  piecewise  linear  representation 
of  the  nonlinear  functional  on  the  goal  arc  is  accomplished.  The  single 
arc  with  nonlinear  functional  between  and  ^ is  replaced  by  a finite 
number  of  arcs  with  linear  functionals  on  each. 


—For  further  detailed  development  of  the  underlying  theory  see  Charnes 
and  Cooper  [ 4 ] Chapter  XVII. 
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FIGURE  2 

GOAL-ARC 

AS 

MULTIPLE  ARCS 


t 


-11- 


Numerical  Illustration 


In  order  to  make  the  preceding  development  more  concrete,  we 
will  now  consider  a numerical  example.  The  problem  that  we  will 
consider  Is  the  problem  considered  by  CharneB,  Cooper,  Lewis  and 
Niehaus.— ^ 

Let  there  be  two  categories  of  personnel  a » 1,  2 (e.g., 
female  and  male)  and  three  time  periods,  t - 0,  1,  2.  For  Job 
categories  we  shall  use  the  following: 


±i_J. 

0 

1 

2 

3 

4 


Description 
Outside  Source 
Clerical 
Technical 
Administrative 
Natural  Attrition 


Abbreviation 

0 

C 

T 

A 

N 


Figure  3 provides  targeted  workforce  goals  a^t)  where 
1 ■ 1,  2,  3 for  the  associated  job  category  in  each  of  the  periods 
t ■ 0,  1,  2.  Figure  4 provides  a matrix  of  transition  probabilities 
which  is  assumed  to  be  applicable  over  these  periods.  Recall  that  N 
refers  to  natural  attrition  so  that,  e.g.,  there  Is  0.26  probability 
that  clerical  personnel  will  leave  the  organization  in  going  from 
one  period  to  another. 


1 


^See  [5] . 


inoum 

1 

2 

c 

B 

700 

650 

T 

875 

450 

400 

A 

225 

200 

200 

FIGURE  3 

TARGETED  WORKFORCE  GOALS,  a^t). 


FIGURE  4 

EXAMPLE  MARKOFF  MATRIX 


-13- 


In  Figure  5 Che  actual  proportions  of  personnel  In  each 
Job  category  for  the  Initial  time  period  and  the  desired  proportion 
of  personnel  In  each  job  category  for  future  time  periods  are  given. 

The  actual  proportions  are  obtained  from  the  "on  board"  starting 
population.  The  desired  proportions  represent  policy  statements 
concerning  the  desired  mix  of  personnel  for  the  future. 

Figure  6 provides  the  desired  number  of  personnel  of  type 
a«l  (female)  for  each  job  category  In  each  period.  These  values  are 
obtained  from  Figure  3 In  the  following  manner.  Let  b^t)  ■ 
where  ^u^  is  the  smallest  Integer  not  less  than  u.  Thus,  e.g..  In  Figure  6 
525  ■ .75  x 700  in  the  row  for  C where  it  intersects  the  column 
captioned  "1"  is  obtained  from  the  data  of  Figures  3 and  5. 


Actual 

Propor- 

tions 


Female 

ON 

00 

.20 

.40 

Male 

.11 

o 

00 

.60 

Female 

.75 

.35 

.45 

Male 

.25 

.65 

.55 

Desired 

Propor- 

tions 


FIGURE  5 

EXAMPLE  OF  PERSONNEL  - JOB  PROPORTIONS,  p“ 


FIGURE  6 

TARGETED  FEMALE  WORKFORCE  GOALS,  b^t) 
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In  Figures  7 and  8 the  "artifact  goals"  are  given  for  each 
, of  the  two  periods  as  Indicated  In  the  titles  of  these  Figures.  The 

s "artifact  goals"  are  defined  by  g“j(t)  ■ ^P*a^(t-1)  where 

Is  the  ltj^  element  of  the  Markoff  matrix  M.  In  this  example  we  are 
confining  our  attention  to  a-1  and  so  we  can  let  8^(0  “ g^Ct) 
without  ambiguity. 

Similarly  let  x^Ct)  equal  the  number  of  females  (a«l) 
transferred  from  job  category  1 to  Job  category  j In  period  t and  let 
yj(t)  represent  the  total  number  of  females  In  job  category  j in  period 
t.  In  this  model  the  y^(t)  and  the  x^(t)  are  to  conform  "as  close  as 
possible"  to  the  targeted  workforce  goals  and  the  "artifact  goals" 
respectively. 
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Reduction  to  Network  Format 

We  now  formulate  this  as  a network  problem.  This  Is  shown 
graphically  In  Figure  9.  In  this  example  J*,  j"(t)  ■ K~  and  V^t) 

■ where  K Is  an  abbreviation  for  "Job  category."  Here,  of 

course,  K takes  on  the  values  A,  C,  T,  N,  0.  As  already  noted,  the 
symbol  ^ on  an  arc  indicates  that  it  is  a "goal  arc."  Upper  and 
lower  bounds  for  the  flow  on  the  "valve"  arcs  are  set,  respectively, 
at  the  projected  manpower  requirements  plus  ten  per  cent  of  the 
requirements  and  minus  ten  per  cent  of  the  requirements. 
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In  this  example  we  will  employ  only  two  pieces  in  our 
piecewise  linear  goal  functional,  i.e.,  k - 2.  Hence  the  decomposition 
on  a "goal  arc"  is  performed  as  described  earlier  with  k = 2.  We  will 
now  examine  the  decomposition  of  "goal  arcs"  in  this  example. 


Consider  any  goal  arc  in  Figure  9 between  a and  a Kt. 

We  replace  this  arc  in  Figure  10  with  two  arcs,  say,  G^(t),  where 

Jq 

k * 1 or  2.  Let  x. . (t)  denote  the  corresponding  flows.  These  flows 

1 2 

are  bounded  as  follows:  0 x(j(t)  < g^.(t)  and  0 < x^(t)  < °°. 

k k 

Let  c denote  the  functional  coefficient  on  G^(t).  We  assume  that 

12  2 
c < c . In  an  optimal  solution  there  wLLi  be  no  flow  on  G^j(t)  until 

the  flow  on  G ^ j ( t)  has  reached  g^(t). 

Now  consider  a goal  arc  between  nodes  K+  and  K , , . As  above, 

t t+1 

12  k 

we  replace  this  arc  with  two  arcs,  G^(t)  and  G^(t).  Let  y^(t)  denote 
|( 

the  flow  on  G . ( t ) . The  flows  on  the  two  arcs  are  bounded  as  follows: 

1 2 k 

'<  ■/  i t > b^(t)  and  0 < y^(t)  < *•<.  Let  d denote  the  functional 

k 12 

coefficient  for  the  flow  on  0j(t).  We  assume  that  d < d . 

Proceeding  in  this  manner  the  problem  is  represented  as  a 
network  with  the  "goal  arcs"  decomposed  as  In  Figure  10. 


i 


NETWORK  WITH  GOAL-ARCS  DECOMPOSED 
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Slnce  the  objective  function  is  to  be  minimized,  a high  positive 
value  for  the  functional  coefficient  on  an  arc  tends  to  make  the 
resistance  to  flow  on  that  arc  high.  In  our  penalty  system  the 
following  priorities  are  established:  Meeting  the  goal  of  a certain 
number  of  female  personnel  for  each  job  category  in  each  time  period 
la  given  the  higest  priority.  Firing  is  highly  discouraged.  Flexible 
movement  has  the  second  highest  priority.  The  penalty  on  exceeding 
manpower  requirements  is  greater  than  any  other  penalty  except  the 
penalty  on  firing.  The  penalty  for  hiring  in  the  first  period  is 
greater  than  the  penalty  for  hiring  in  the  second  period.  The  penalty 
on  hiring  is  less  than  the  penalty  on  exceeding  manpower  but  greater 
than  the  penalty  on  flexible  movement.  The  penalty  on  firing  is  set 
at  an  order  of  magnitude  larger  than  the  sum  of  all  other  weights. 

The  values  for  the  functional  coefficients  on  the  arcs  (with 
relevant  interpretations)  are  given  as  follows: 

H = hiring  penalty  = 5; 

P - penalty  on  flexible  movement  = 2; 

R - firing  penalty  = 1,000; 

G - penalty  on  expected  movement  = -1; 

Q » penalty  on  meeting  manpower  requirements  = -6; 

F - penalty  on  exceeding  manpower  requirements  • 10. 

The  solution  is  summarized  in  four  tables  as  follows:  The 

projected  personnel  transfers  for  periods  1 and  2 are  given  in  Tables  1 
and  2,  respectively.  The  424  under  "Normal  t Flexible"  in  row  1 of 


I 


i 


I 

i 

I 

! 


I 

I 

I 


Table  1 represents  the  planned  retention  of  females  In  the  clerical 
Job  category  in  the  first  time  period.  It  is  composed  of  420  females  via 
normal  retention  plus  4 more  as  a part  of  an  optimum  managerial  plan 
to  alter  the  present  composition  of  the  organization.  The  total 
of  525  females  at  the  bottom  of  this  column  is  to  be  obtained  by 
recruiting  an  additional  101  females  from  outside  the  organization. 

Table  2 is  similarly  interpreted  for  the  second  time  period. 


I 

I 


PROJECTED  PERSONNEL  TRANSFERS 
FEMALE 
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Table  3 compares  workforce  requirements  and  the  optimal 
distribution  from  the  model  — e.g.,  targeted  workforce  goals  and  optimal 
"aboards."  The  discrepancies  between  the  two  are  given  in  the  last 
column  of  Table  3.  All  discrepancies  are  at  zero  value  which  means 
that  the  optimum  program  achieves  all  of  the  indicated  targets. 


TABLE  3 
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Table  4 is  a summary  of  the  personnel  actions  projected  by 
the  optimum  plan.  For  example,  420  normal  transfers  plus  4 additional 
(flexible)  transfers  and  101  hires  are  projected  for  the  clerical 
category  in  Period  1 and  368  normal  transfers,  2 additional  (flexible) 
transfers  and  118  hires  in  period  2. 

Summary  and  Conclusion 

This  concludes  the  present  paper,  but  the  above  developments 
are  a continuation  of  research  in  a series  dealing  with  modeling  for 
EEO  planning.  The  first  in  this  series  of  models  was  the  FEEO  model 
which  provides  for  EEO  planning  at  the  macro-level.  See  [6].  The 
next  in  the  series  was  the  MEEO  model.  Also  called  the  "Coherence  Model," 
the  MEEO  model  was  developed  to  provide  for  EEO  planning  at  the  micro- 
level, e.g.,  at  the  activity  level,  which  would  be  "coherent  with"  the 
FEEO  model.  For  further  discusrion  see  [5]. 

The  model  developed  above  is  an  extension  and  reformulation  of  the 
MEEO  model.  As  such  we  have  a continuing  evolution  in  a modeling  strategy. 
The  problem  which  was  originally  formulated  as  a capacitated 

distribution  problem  with  "artifact  goals"  is  now  reformulated  as  a network 
problem  with  goal  arcs.  Thus  we  have  alternate  models  for  this  same 
class  of  problems.  The  development  portrayed  in  this  paper  was  undertaken 
to  take  advantage  of  large  capacity,  fast  and  highly  efficient  network 


codes  such  as  PNET.  See  [7]  and  [8].—^  Also  In  the  model  developed  above 
the  transhipment  characteristics  provide  much  greater  convenience, 
simplicity  and  flexibility  in  representing  desired  personnel  flows.  In 
the  MEEO  model  the  dyadic  character  required  special  devices  and  redundant 
representation. 

The  Coal-Arc  model  of  this  paper  currently  handles  the 
ethnosexual  categories  one  at  a time.  This  is  done  via  the 
proportionate  reduction  devices  described  in  the  above  paper.  The 

( 

next  step  in  this  ongoing  research  should  develop  a method  for 
handling  all  of  the  ethnosexual  categories  simultaneously.  This  and 
other  parts  of  this  work  in  EEO  modeling  will  be  reported  in 
subsequent  papers  of  this  series. 


—See  also  [1]  and  [2]  and  [9]. 
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